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We provide the most general forms of covariant and normalized time operators and their prob- 
ability densities, with applications to quantum clocks, the time of arrival, and Lyapunov quantum 
operators. Examples are discussed of the profusion of possible operators and their physical meaning. 
Criteria to define unique, optimal operators for specific cases are given. 
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I. INTRODUCTION 

Quantum theory provides, for a given state prepara- 
tion, expectation values and distributions for a number 
of observables whose operators have been identified by 
a combination of heuristic arguments (e.g. quantization 
rules), and consistency arguments. Their relevance and 
validity is eventually put to the test or motivated by 
experiments. Time observables, i.e., random variables 
such as the arrival times of particles at a detector for 
a given state preparation, have been more problematic 
than other observables like "energy", "momentum" or 
"position" evaluated at a fixed instant. In fact almost a 
century after the creation of the basic quantum formal- 
ism, the theoretical framework to deal with time observ- 
ables which have a relatively straightforward operational 
definition in the laboratory, is still being debated. Re- 
views of several aspects of the difficulties and efforts to 
formalize time in quantum mechanics may be found in 
two recent books [l|, [2[ . Some of these difficulties may be 
traced back to a lack of a general framework to generate 
and define "time operators" . An important point, fre- 
quently overlooked, is that, for a given system, there is 
no single time operator. There are infinitely many time 
operators corresponding to different observables and ap- 
paratus. "Canonical time operators" have been defined 
0,3, but, as we shall stress, the definition of "canoni- 
cal" is basis dependent, even without energy degeneracy. 
Thus, further analysis is necessary to set ideal operators 
and possibly uniqueness in some cases by imposing the 
physical conditions to be satisfied (e.g. symmetries) or 
optimal properties, such as a minimal variance. 

Time operators can be classified into two main groups 
on physical grounds, depending on their association with 
time durations or time instants. An example of a dura- 
tion is the dwell time of a particle in a region of space. 
The corresponding operator commutes with the Hamil- 
tonian since the duration of a future process does not 
depend on the instant that we choose to predict it |5|- 
Instead, the other group of time observables are shifted 
by the same amount as the preparation time, either for- 
ward (clocks) [6[ or backward (event times recorded with 
a stopwatch, the simplest case being the time of arrival) , 
and are conjugate to the Hamiltonian. We shall set here 



a framework for these "covariant" observables [3[ asso- 
ciated with instants and analyze their multiplicity and 
physical properties. Applications are discussed for quan- 
tum clocks and the time of arrival. The relation to Lya- 
punov operators is also spelled out. 

The plan of the paper is as follows. After introducing 
the main concepts and notation in Sec. |lll in Section 
mil the most general form of a covariant time operator is 
determined for a Hamiltonian with only continuous, pos- 
sibly degenerate, eigenvalues. In Section IIVI it is shown 
that, for a time-reversal invariant Hamiltonian, one ar- 
rives at a unique and natural form of time operator by 
imposing time reversal covariance, invariance under ad- 
ditional symmetries and minimality of the variance. In 
Section |V] the results are applied to arrival times for a 
particle moving on a half-line and a connection with the 
delay time of Smith Q is established. In Section IVll the 
results are applied to Lyapunov operators which were 
considered in Ref. Q- It is shown that the expression 
given there is a special case, and the general form as 
well as possible uniqueness conditions are presented. In 
particular it is shown that for a time reversal invariant 
Hamiltonian there is no time reversal invariant Lyapunov 
operator. This is of interest because it has been argued 
that, in order to characterize a quantum system as irre- 
versible and an arrow of time if the Hamiltonian is time 
reversal invariant and if one uses a formulation in terms of 
Lyapunov operators, a Lyapunov operator or functional 
should be time reversal invariant |9| . 



II. COVARIANCE OF TIME OPERATORS. 
NOTATION. 

We differentiate between clock time operators and 
event time operators. The former, denoted by T, can 
be associated with a quantum clock which measures the 
progressing parametric time, while the latter, denoted by 
T , describe the time of an event, for example, the in- 
stant of time a particle is found to arrive at a particular 
position. This and the following two sections are mostly 
devoted to clock observables, although the formal results 
are analogous for event times. In an ordinary clock the 
dial position is the observable which tells us what time it 



is. In a quantum clock the dial "position" is probabilistic 
but its average should follow faithfully the advancement 
of parametric time. Wc would like as well to minimize 
the variance and estimate the time as accurately as pos- 
sible with a finite number of measurements. We will not 
investigate here specific operational realizations, see a 
review in |6|, but instead idealized operators and their 
properties. 



A. Clock time operators 

For a given state \ip), let the probability of finding 
the measured time in the interval (— oo,r) be given by 
the expectation with Jt/i) of an operator Fr- Note that 
< Ft- < 1 so that Ft is sclfadjoint and bounded. [For a 
momentum measurement the analogous operator would 
be J_ dp'\p'){p'\ for finding the momentum in (— oo,j)). 

Here F-r can have a more general form and in general one 
deals with a positive-operator valued measure.] Then 



Ii{T-ij) = —{i,\FM) 



(1) 



is the corresponding temporal probability density, nor- 
malized as j dTli{T]ip) = 1. We define the probability 
density operator IIt- by 



normalized as 



n^ " Tr^- 



drUr = 1. 



(2) 



(3) 



The mean value of observed time can be written as 



dTTll{T;^) = {iPl dTTllr\^) = {^\f\ij), (4) 



where 



T : 



dr rllr 



is called the time operator associated with IIt-. The sec- 
ond moment, if it exists, is given by 



dr T^ n(T; V) = (^1 / dr r^ H^ |^) 



(6) 



and similarly for higher moments. It may happen that 
this is not equal to {ip\T^\ip) . 

A clock time operator is called covariant with respect 
to ordinary (parametric) time if for the states \ip) = jV'o) 
and \ipt) the probabilities of finding the measured time in 
the respective intervals (—00, r) and (—00, r+t) coincide, 
i.e. if 



This implies, choosing t = — r. 

Ho = ^[^,^0] (9) 

n^ = e-'^^/^tloe'"^^". (10) 



Note that (?/'|IIr|'0) is non-ncgativc because {tp\FT\tp) is 
non-decreasing. By a change of variable in Eq. ([S]) one 
obtains 



iHt/h rp —iHt/h rp I . 



(11) 



From this it follows by differentiation that H and f sat- 
isfy the canonical commutation relation 

[f,H]^ih (12) 

when sandwiched between (normalizable) vectors from 
the domain of H. 

Note that Hq and IIt- arc in general not operators on 
the Hilbert space but only bilinear forms evaluated be- 
tween normalizable vectors from the domain of H. An 
expression like {E\IIo\E') has to be understood as a dis- 
tribution. Since the diagonal E = E' has measure it 
is no contradiction that Eq. ^ gives on the diagonal 
while the following example gives (27rfi,)~^. 
Example: For a Hamiltonian H with non-degenerate 
continuous eigenvalues E and eigenvectors \E) with 
{E\E') = 5{E - E') we put 

{E\Ilo\E')^ ^ 



so that in this case 
1 

27rn 
1 

27r^ 



Ho = 



n. = 



dE dE'\E){E'\, 

dE dE'e-'^^~^'^^^''\E){E'\ 



(13) 

(14) 
(15) 



The normalization condition of Eq. ^ is easily checked. 
The corresponding clock time operator —ihdE which re- 
sults from Eq. ([5]), has been considered the "canonical 
time operator in the energy representation" 
but note that | E) is unique only up to a phase 



(5) and taking \E)^ 



^ piviE) 




\E) instead of \E) leads, for dif- 



ferent If, to multiple "energy representations" , even for 
a system without any degeneracy. In the new basis the 
"canonical operator" will be shifted by 



h / dE ip'{E)\E){E\. 



(iPolFMo) = {A\Fr+t\A) 



(7) 



(16) 



Moreover, the mean-square deviation AT^ for a given 
state depends on (p{E) in such a way that there is no 
choice of (p{E) which would make AT minimal for all 
states, as shown in Appendix El Therefore, in this case a 
minimality condition imposed on AT cannot be fulfilled 
and does not lead to a unique natural choice of time 
operator without further additional restrictions. There 
must be additional physical criteria to choose, and in 
fact several of them may be physically significant. This 
will be exemplified below, see Sect. |Vl 



B. Arrival time operators 

"Time-of-event" , and in particular timc-of- arrival op- 
erators and probability densities are similar to clock oper- 
ators (for reviews of this concept see [ij,[ij|). Physically, 
we expect that a free particle in one dimension will ar- 
rive with certainty at a given detection point (including 
negative times and ignoring the case of zero momentum 
which is of measure zero for an arbitrary physical wave- 
packet). Similarly a free particle in three dimensions will 
arrive at an infinitely extended plane. Also, a particle on 
a half-line with reflecting boundary conditions and with- 
out additional potential, is expected to arrive once at the 
boundary and, at least on classical grounds, twice at any 
other point. In the latter case it is meaningful to consider 
the first arrival at a given point because this should be 
in principle observable. In all these cases the total ar- 
rival probability resp. first-arrival probability is 1. The 
corresponding arrival time operators are denoted by T 
and n^, respectively, and when compared to clock oper- 
ators their formal properties are identical up to a change 
of sign, e.g. in the conjugacy relations or the formu- 
lation of covariance [T^- This means that, in contrast 
to clock times, if the particle's state is shifted in time 
by to, it should arrive a time to earlier, and the tempo- 
ral probability density should be shifted by to to earlier 
times. These are, in other words, waiting times until an 
event occurs, which depend on the time when we set the 
stopwatch to zero, and decrease if we reset it at a later 
instant. Thus the analog of the cumulative probability 
operator in Eq. ([7]) must now satisfy 



(V'ol^.^l^o) = {^MFr^-tHt) 



F' 



iHr/h fp/ 



-iHr/h 



(17) 



With n^ = dF^/dt and fl^ = \[H, F^] we have 




f^ = f dtte'^'/^ti^e- 


iHt/h 

1 


(18) 


nf = e^^*/''fI^e-'^*/^ 




(19) 


(^tolT^l^to) = (^o|T^|^o)-<o, 




(20) 


(V-toinfl^to) = (^oin^+tJV'o). 







In addition, the operator should incorporate the loca- 
tion where the arrivals are observed. For free particles 
coming in from one side and arrivals at a plane this was 
achieved in Ref. [16[ by postulating invariance of the 
probability density under a combination of space reflec- 
tion and time reversal. It is evident that these properties 
still do not specify the operator uniquely. For physical 
reasons one will also demand for an optimal arrival-time 
observable that the arrival-time probability density has 
minimal variance, analogous to the postulate in Ref. [16[ 
for free particles in three-dimensional space. This means 
that no other arrival-time observable can be measured 
more precisely. 



III. THE GENERAL FORM OF COVARIANT 
TIME OPERATORS 

We begin with covariant clock time operators associ- 
ated with a given Hamiltonian H . For simplicity, we first 
consider the case when H has only non-degenerate con- 
tinuous eigenvalues E, with generalized eigenvector \E) 
and normalization 

{E\E') =6{E-E'). 

We will determine the most general form of IIo which, 
through Eqs. ([T]-[TU]), leads to a covariant probability 
density operator and corresponding time operator. 
The simple example in Eq. (|14p can be generalized to 



1 



no^TTT dEdE'b{E)\E){E'\b{E') 
Znh. J 

and. more generally, it will be shown that 

no - ^E/^^^^'^'(^)I^)(^'I^^(^' (21) 

i 

f = J^J2 fdtt fdEdE' e-(^-^')*/'"' 

X b,{E)\E){E'\b;(E^ (22) 

is the most general form of IIq and T, where the functions 
hi [E) have to satisfy certain properties in order that the 
total probability is 1 and that the second moment in 
Eq. ([6]) is finite. Indeed, for given state |-0), the total 
temporal probability is, with V'(-£') = {E\i'), 



+ 00 



dt(V|e-*^*/'^noe*^*/''|V) 



(23) 



E 



dt 

2TTh 



dEe~'^''^TP{E)b,{E) 

"E j dEdE' 5{E - E')iJ{W)h[E)b~(W)i){E') 

= Y, I dEi;{E)Y,h,{E)h{E)^{E). 

This equals 1 for every state jf/;) if and only if 

Y,h{E)h;{E) = ^- (24) 

i 

Similarly, 

(V'iflV) = ^dEi,{E)-^^'{E) 

+ JdE\^{E)\^^Y.^,iE)bm. (25) 

Note that ^ bi b[ is purely imaginary, from Eq. ^IQ, and 
thus vanishes if 6,- is real. 



The second moment is 

dt 

2^ 



(26) 



= ^^E / ^ / dEdEe-'''"''^{E)h{E) 

i 

= h^Y.h^^E [m) h{E) ) dE ih~(E)i;{E)^ 
+2 Re ^6^6^(i?)^K^V'(£;) I 



by Eq. ([24|) . This is finite if and only if the contribution 
from the first and second term arc finite, and for the lat- 
ter to hold for all infinitely differentiablc functions ipi^) 
vanishing outside a finite interval (i.e. with compact sup- 
port in E) one must have 

\^|6^(£')P integrablc over any finite interval. (27) 



Eq. (PT|) gives the most general form of XIq leading to 
a covariant time operator when the functions bi satisfy 
Eqs. (p4)) . and the second moment is finite for states with 
{E\il}) of compact support if and only if Eq. (P7|) holds. 
For a given Ho one can construct the functions bi as 
follows. One chooses a maximal set {|(7i)} of vectors sat- 
isfying 



(.9i|no|gi) = 5ij/2Tih. 



(28) 



Such a maximal set is easily constructed by the standard 
Schmidt orthogonalization procedure. Then a possible 
set {bi\ is given by 



b,{E)^2nrl{E\t\o\g^ 



(29) 



Eq. (PT|) is then a realization of the given XIq. Math- 
ematical details, in particular regularity properties, will 
be presented elsewhere [l^. It should be noted that the 
functions bi in the decomposition of Ho in Eq. (PT|) are 
not unique. 

For the case of degenerate eigenvalues of H we first 
consider the case where the degeneracy is indexed by a 
discrete number and such that 



{E,a\E',a') ^^^ Saa'S{E - E'). 



(30) 



For simplicity we assume the same degeneracy for each 
E. Then Eqs. ([21] -[211) generalize as 



""-^E 



dEdE' 



Y,HE,a)\E,a){E',a'\b,{E',a'), 



(31) 



Y,h{E,a)b,{E,a')^Sa 



(32) 



second moment ~ h^ dE Ids /^ bi{E,( 



.)ijiE, 



(33) 



y \bi{E, a)p integrable over any finite interval (34) 

i 

for each a, where ip{E,a) = {E,a\ip) and bi{E^a) = 
27r?i(£',a|IIo|(7i). Again Eq. ([3T|) gives the most general 
form of Hq leading to a covariant time operator through 
Eqs. (fniTU)). The case of continuous degeneracy parame- 
ter can be reduced to the discrete case. 

These results carry over in a corresponding way to ar- 
rival times with normalized probability densities. 



IV. UNIQUENESS OF TIME OPERATOR: TIME 

REVERSAL, SYMMETRIES AND MINIMAL 

VARIANCE 

As seen in the previous section, there are many co- 
variant clock time operators. For uniqueness additional, 
physically motivated conditions are needed. Requiring 
minimal variance by itself does not make T unique, not 
even in the case of non-degenerate spectrum of H, since 
in general it may not be possible to fulfill this require- 
ment simultaneously for all states with second moment 
unless, in addition, one restricts the set of functions bi 
by symmetry requirements, as we shall now discuss. 

The time reversal operator, here denoted by Q, is an 
anti-unitary operator. If the dynamics is time reversal 
invariant, it is natural to demand 



ere = -T, 



(35) 



and similarly for the probability density. By Eq. ([9]) this 
implies 



enoe = no. 



(36) 



It will now be shown for the non-degenerate eigenvalue 
case that time reversal invariancc of the Hamiltonian H 
and of Ho, and minimal AT together imply uniqueness 
of T and lit • For each eigenvalue E oi H one can choose 
a Q invariant eigenvector, denoted by \Eq), 



eiEf.) = \Ep, 



(37) 



This means a specific choice of phase factor and a real 
function in position space. Eq. p6p implies Hq = 
l/2(IIo + 9IIo0), and the general form of Ho in Eq. 
(PT|) then implies that bi{E) can be chosen real. Then, 
from Eqs. (^5)) and ^Q, one finds 



(^|f|^) 



dEME)-ij'(E), 
I 



(38) 



second moment = h"^ J dE |^'(£')p 

+h'E^IdE\^P\'mE)\\ (39) 

Hence AT minimal means in this case that the second 
moment is minimal, and the latter holds if and only if 
briE)^0,i.e. 

h{E) = Q, 5] c2 = 1, 

by Eq. (f24| . Inserting this into Eq. (f2T|) one sees that 
the functions bi can be replaced by the single function 
h{E) = 1. Thus one obtains 

1 



(40) 



with time reflection invariant \Eq). The (non- 
orthogonal) eigenfunctions, |r), of T with eigenvalue r 
are given by 



n^^^^j dEdE'\Ee){E'Ql 

n* = ^ J dEdE' e-^(^-^'^*/^\Ee){E'^l 



T = dtn 



dEe-''^^'^\Ep 



and T can be written as 

/oo 
drr|r)(r|. 
-OO 



(41) 



(42) 



Therefore uniqueness holds in the non-degenerate case if 
time-reversal invariance and minimal AT are demanded. 
In the degenerate eigenvalue case this is no longer true 
and one needs additional conditions to obtain unique- 
ness, as discussed elsewhere [I4]. Here we simply state 
some results. With a reflection invariant potential in one 
dimension, the clock time operator becomes unique and 
can be explicitly determined if, in addition to covariance 
under time reversal and minimal variance, one also de- 
mands invariance under space reflection. With a rotation 
invariant potential in three dimension, the time operator 
becomes unique and can be explicitly determined if, in 
addition to covariance under time reversal and minimal 
variance, one also demands invariance under rotations 
and reflection xi — >■ —xi. Analogous results hold for ar- 
rival time operators. In particular, a generalization of 
the result of Ref . [l6[ is obtained [13| • 



V. APPLICATION TO ARRIVAL TIMES 

Evidently the techniques of the previous sections can 
be applied in a completely analogous way to the study 
of arrival-time operators. To illustrate this we consider 
in the following the motion of a particle on the half-line 
a; > 0, without additional potential, and study its arrival 
times at the origin and at an arbitrary point. 



In the classical case an incoming free particle of en- 
ergy E is reflected at the origin and then travels back 
to infinity. Hence, for each point a 7^ 0, there is a first 
and second time of arrival which we denote by i° and 
^2- For the time reversed trajectory the first arrival at 
a is at time ig ^ = —f^ and the second arrival at time 
ig 2 = ~tl^ s-s is easily calculated. For the origin, a = 0, 
there is only one arrival and 






-t°. 



(43) 



The corresponding arrival-time operator for arrivals at 

f 



the origin is denoted by Tf. It is natural to demand the 



analogous relation to Eq. (|43p . i.e. 



ef/e 



rfiA 



(44) 



and time reversal invariance of Hf q' 



where Yl4^ is the 



associated probability density operator. 

If a 7^ a classical free particle on the positive half-line, 
coming in from infinity with velocity |w|, arrives first at 
time ti at the point a, and then at time t*^ at the origin. 



t^ = t°-a/\v\. 



(45) 



If T{^ denotes the corresponding time operator for the 
first arrival at a we may demand 



Tf 



Tf-a/\v\ 



(46) 



where \v\ = \/2H/m is the velocity operator. 



A. Free particle on a half-line 

We first consider arrivals at the origin for free motion 
on the half-line x > 0, with reflecting boundary condi- 
tions at a; := 0. The eigenfunctions can be labeled by 
the energy E = k'^h?/{2m). Real, and thus 8 invariant, 
eigenfunctions for energy E which vanish at the origin 
are 



(rlEf 



27rfc 



(e 



'ikr 



(47) 



where the subscript / in \Ef) refers to the free Hamilto- 
nian and where we have written r to indicate r = x > 0. 
These eigenfunctions are normalized as {Ef\E'f) = 5{E — 
E') on the half-line. 

For the probability density operator for arrivals at the 
origin invariance under time reversal means 



&tV%Q 






(48) 



By the results of the last section, the operators H^^ and 
T^ become unique if invariance under time reversal holds 
and minimal variance is assumed. From Eq. (|40p one 



obtains, with a change t -^ —t and replacing \Eq) by 



n 



/'* = ^JdEdE'e^(^-^'^'/'^\Ef){E'f\, (49) 
ff = jdttlil. 



This arrival time operator is just the negative of the clock 
time operator of Eq. (gD]), with Eqs. (gT]) and (g^]) hold- 
ing correspondingly. 

Note that the vanishing of the wave function at r = 
is not an obstacle to define these operators in a physically 
meaningful manner. A similar situation is found for anti- 
symmetrical wavefunctions on the full line. It was shown 
in [13I that the ideal timc-of-arrival distribution follows 
in a limiting process from an operational measurement 
model that considers explicitly a weak and narrow detec- 
tor. 

We now turn to first arrivals at a 7^ 0. Using Eq. [T2| 
a simple calculation shows that 



\v\lh rp 



f 



'1*1/1 



T- 



/ 



(50) 



Since the right-hand side equals f^, by Eq. p6)) . this 
implies an analogous relation for the probability density 
operator, 11^ j, for T^, 



(51) 



TjA iam\v\/hY\A —iarrt\v\/h 

Using Eq. (gHl) this can be written as 

(52) 
which explicitly gives the temporal probability density 
operator for the first arrival at the point a of a free par- 
ticle on the positive half-line. For a — ?► one recovers 
Eq. (USD. 



The M0llcr operators Ct± satisfy 

/•CXD 

n± = lim e'"*^^'e-'^f*^'' = / dE\E±){Ef\, 

\E±) = n±\Ef). (54) 

The freely moving asymptotic states I'lpin) and \^pout) are 
mapped by fl± to the actual state \ip)^ 



IV) = r!±|V',„) 

out 
\i>out) = S\lpin) 



(55) 



where S — Ol_^+ is the S operator. Note that, by 
Eq. dSSD, 



(56) 



S= dE\Ef)e'^\Efl 

Jo 

so that e^^^ is the eigenvalue of S for the state \Ef) 
It is convenient to introduce also the operator 



ne= / dE\Ee){Ef\ 
Jo 



and define operators T^q by 



T. 



±,e 



^±,eTf ^j-Q 



(57) 



(58) 



= fdttfdE dE' e*(^-^')*/'^|^±,e) {E'^ 



±,ei- 



The last line shows that —T_^q are possible clock time 
operators for the particle in the potential. Since the 
states \E±q) differ only by a phase, the same calcula- 
tion that leads to Eq. dTS]) gives 

f^=f,^ThfdE^\Ee){Ee\. (59) 



B. Asymptotic states and Smith's delay time 



We now apply the free-particle result in Eq. (|49)) to 
the asymptotic states of a particle in a potential on the 
half-line whose Hamiltonian has no bound states and to 
which scattering theory applies. Although for fixed E 
the eigenstate is unique up to a phase, there are phys- 
ically relevant eigenstates \E±) which correspond to an 
incoming (-1-) and outgoing (-) plane wave, respectively, 
as well as the invariant state, denoted by |£'e)- Their 
relation and asymptotics are |£'-) = B|i?+) and, with 
the scattering phase shift S = 6{E), 



{r\E+) 

{r\E-) 
{r\Ee) 



1 /2m i _jj,^ 



J2i5 Akr\ 



ME. 



e-'^'Hr\E. 



+ h 



e-'Hr\E, 



(53) 



From Eq. ([55|) it follows that the expectation values of 
T+ , T_ and Tq may be interpreted in terms of the 
asymptotic states and the free-motion arrival-time oper- 
ator ff, 

(■(/'|T'_/^_,elV') = {'^l^in.oiit,%o\Tf\i'in,out.io), (60) 

where the freely moving state |Vio) is defined by 

|V„) = 5i/2|^„), (61) 

and can be considered as an interpolation between \ipin) 
and IV'out) = S\ilJin)- With Eq. ((57)) one can write 



\i^^o)=^lH). 



(62) 



Taking expectation values of Eq. ([5^ with \il}) and us- 
ing Eqs. (jnni) and ([M]) . together with the fact that 



\E±.e){E±,Q\ all coincide since the phases drop out, 
yields 



dE 



(63) 



One sees from this that the mean arrival time for the 
interpolating state IV'jo) lies between those of the ingoing 
and outgoing wave. From Eq. (|63p . 



dE^\{Ef\i^,.^ 



{1i^out\fJ^\^l^out)-{'>i^^u\ff'^\'<P^n) = iH 



(64) 
The right-hand side of the last equation is the scattering 
time delay of Smith Q and it shows that the time for the 
outgoing wave is shifted with respect to the time for the 
ingoing wave by the scattering time delay. An example 
is shown in Figs. 1 and 2. 

Time-reversal: The behavior of T^ with respect to 
time-reversal is determined by acting with the anti-linear 
operator 9, 



erre^ ~T4 



(65) 



whereas Tq simply changes sign. The operators T^ do 

not simply change sign under time-reversal as Tq does, 
but their behavior in Eq. (|65p (changing the sign and 
exchanging the operators) is perfectly physical: the time 
reversal of a trajectory which moves towards the origin is 
a trajectory in the same location but moving away from 
the origin. If the original incoming trajectory requires 
a certain time t to arrive at the origin (with free mo- 
tion), the reversed trajectory is outgoing, and departed 
from the origin at — r. These operators provide in sum- 
mary information of the free-motion dynamics of incom- 
ing and outgoin g as ymptotes of the state, and scattering 
time delays [13,123 ■ Thus, although the operator Tq is 
unique when one applies the criteria of the previous sec- 
tion it does not supersede T^ since it does not describe 
the same physics, and all three operators have their own 
legitimacy. 





FIG. 1: (Color online) Probability densities before (i = 
0) and after the collision (i = 190) with a delta bar- 
rier. Dimensionless units with m = h = 1. The ini- 
tial wave packet is '(/'(fc) = A''[l — exp(— /3A:^)] exp[— (fc — 
fco) /(4A(,)] exp(— i/sa::o)^(fc), where A'' is the normalization 
constant and 9 (here) the Heaviside step function; initial 
wavenumber ko = -7r/2, Afc = 0.045, {3 = 1/2; V = 
205(2: — 20); initial center of the wave packet xo = 180. The 
delta potential is rather opaque so the the outgoing packet is 
advanced with respect to the incoming state. 
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FIG. 2: (Color online) Time of arrival distributions for ar- 
rivals at X- = corresponding to the previous figure. 



VI. APPLICATION TO LYAPUNOV 
OPERATORS IN QUANTUM MECHANICS 

In Ref. Q an operator L was called a Lyapunov oper- 
ator if for any normalized \ip) and \^pt) = e~*^*/''|V'), the 
expectation value {'ipt\L\4't) is monotonically decreasing 
to as i —> oo and goes to 1 for t — > — oo. Ref. Q consid- 
ered the case of a Hamiltonian H with purely continuous 
eigenvalues ranging from to infinity and degeneracy pa- 
rameter j. The particular Lyapunov operator suggested 
there can be written as 



L, 



3 



dE 



dE' 



AE,j){E\j\ 



E-E' 



■IE 



(66) 



More generally, one may call a bounded operator L a 
Lyapunov operator if (■(/it|L|V't) is just monotonically de- 
creasing, without specifying limits. However, it will be 
shown below, after Eq. (p7)) . that, without loss of gener- 
ality, one can always assume the above limiting behavior 
from 1 to as t goes from — oo to -l-oo. 

The above notion does not quite correspond to Lya- 
punov functionals used in Ref. |9| to define irreversibility 
and an arrow of time, since time reversal invariance of the 
functional was assumed there in order to have neutrality 
with respect to past and future. It will be shown further 
below that there are no time reversal invariant Lyapunov 
operators if the Hamiltonian is time reversal invariant. 

It is clear that the above properties do not define L 
in Eq. (j66p uniquely. For example, one can introduce 



phases and still get a Lyapunov operator. In this section 
we are going to determine the most general form of L for 
a Hamiltonian H with a purely (absolutely) continuous 
spectrum and give conditions under which it becomes 
unique. It will also be seen that to each L there is an 
associated covariant time operator Tl- 

To show that one can assume the above limit behavior 
we put, for a given general Lyapunov operator L, 



Lt = e-'"'/^U"''^ 



(67) 



so that Lt is monotonically increasing, by the monotonic 
decrease of (■0t|i|V't)- From the boundedncss of L and 
from monotonicity it follows that L±oo exist as opera- 
tor limits in the weak sense, i.e. for expectation values. 
Moreover, i±oo commutes with e^*^*/'', and therefore 
i' = L — £_oo is also a Lyapunov operator, with L'^ > 0. 
Then L" = L'^^'^L'L'^^'^ is a Lyapunov operator satis- 
fying i"oo = and L'^ = 1 so that {tljt\L"\^pt) is mono- 
tonically decreasing from 1 to 0, proving the above claim. 
To determine the general form of L with such a limit 
behavior for t — > ±c», we note that by monotonicity 

flf = ^Lt = e''"'/''^[H, Ly"'"-' > 0, (68) 

at n 

i.e. expectation values of Lt are non-negative for all t, in 
particular 



ii^ = —[H,L]>Q 



(69) 



where the commutator is again to be understood in the 
weak sense via matrix elements and where Ilg is in gen- 
eral not an operator but only a bilinear form, as in Eq. 
((9]). From Eq. (|68)) and from L_oo ~ one obtains 



L 



dte-'"*/^tl^e'"'^''. 



From Eq. (|55)) one sees that 

nL{t;ij) = {mf\ij) > 



(70) 



(71) 



is a non-negative density which integrates to 1 for each 
normed state, i.e. it can be regarded as a probability den- 
sity and hence Lt behaves like the cumulative probability 
operator F^ in Eq. ([T|). Therefore, 



fr = / diie-*^*/'*n^e*^*/'* 



(72) 



is an analog of the time operator T in Eq. ([5|) . Alterna- 
tively, 1 — I/_t behaves as the cumulative arrival proba- 
bility operator P/^ in Eq. ([T7]). 

Example: Let Hg given by Eq. (fT4|) . Then, by Eq. 
(ffO)) . L is given by 

L^-^l dt f dE dE'e-'^^-^'^'^^\E){E'\, (73) 



which is readily seen to agree with Ls in Eq. (j66p in the 
case of non-degeneracy. 

For free motion on the half- line, with \E) ~ \Ef) from 
Eq. (|47p , the Lyapunov property of this example simply 
reflects the monotonous accumulation of arrivals at the 
origin since a change of integration variable gives 



(V'tli-ilVt) 



dt'imit'ii^)- (74) 



With a potential on the half-line and taking \E) = \E±) 
of the previous section one obtains the accumulation of 
arrivals of the freely moving packets {ipin) and {ipout), 
and for \E) = \Eq) the corresponding accumulation of 
arrivals for \'>pio)- 

The most general form of L is obtained from the most 
general form of IIq which is given by Eqs. ([31]) and (|32p . 
If IIo is known then L is given by Eq. ([701) , and in this 
way one obtains the most general form of the Lyapunov 
operator L with the above limit behavior for t — > ±cxj. 
Uniqueness of L may be achieved for particular Hamil- 
tonians by demanding, e.g. time reflection invariance of 
Tl, special symmetries and minimal variance AT^, as in 
Sections HV] and |Vl 

We finally show that for a time reversal invariant 
Hamiltonian there is no nontrivial time reversal invariant 
Lyapunov operator. Indeed, ii Q H Q ^ H and L G) = 
L then one obtains, for initial state & IV') = [(OV')), 



{ie^)t\L\ie^)t) = {^-t\L\^-t) 



(75) 



by the anti-unitarity of Q. Now, for increasing t, the 
expression on the left-hand side decreases while the one 
on the right-hand side increases. This is only possible 
if both sides are constant in t. Alternatively, one can 
conclude from Eq. ^ that both 11^ and 9 H^ 9 = -LI^ 
are positive operators, which is only possible if IIq = 0. 
This means that L commutes with H, which also leads 
to the constancy of both sides in Eq. ([75)) . 



VII. DISCUSSION AND OUTLOOK 

We have provided the most general form of covariant, 
normalized time operators. This is important to set a 
flexible framework where physically motivated conditions 
on the observable may be imposed. The application ex- 
amples include clock time operators, time-of-arrival op- 
erators and Lyapunov operators. 

Experimentally, a number of interesting open ques- 
tions remain for quantum clocks and arrival-time mea- 
surements. For example, quantum clocks are basically 
quantum systems with an observable that evolves lin- 
early with time. To evaluate the possibility to compete 
with current atomic clocks |21[, the observable must be 
realized in a specific system. We have described an ideal 
observable (by imposing antisymmetry with respect to 
time reversal and minimal variance) and the analysis of 



the operational realization is now pending. A similar 
analysis for the ideal arrival time-of-arrival distribution 
of Kijowski has been carried out in terms of an opera- 
tional quantum-optical realization with cold atoms (cf. 
Ref. [ij] for a review). Indeed, cold atoms and quan- 
tum optics offer examples of times of events (other than 
arrivals), such as jump times, excitation times, escape 
times, admitting a treatment in terms of covariant ob- 
servables. Modeling and understanding these quantities 
and their statistics may improve our ability to manipu- 
late or optimize dynamical processes. 

On the theory side, an open question is how to adapt 
the proposed framework, possibly in combination with 
previous investigations p^l - l20l [2a - l25l |. to arrival times 
when a particle moves in a potential. 

Finally, we have shown that Lyapunov operators fol- 
low naturally from covariant time observables. Associ- 
ated to time-of-arrival operators, they account for the 
monotonous accumulation of arrivals for freely moving 
asymptotic states from the infinite past independently of 
the state chosen. Note that the "infinite past" here is 
an idealized construct since it must be assumed that the 
wave has been evolving forever, ignoring the fact that in 
practice the state may have been prepared at some spe- 
cific instant. In other words the Lyapunov operator does 
not depend on that preparation instant, and when ap- 
plied to the state it takes into account its idealized (not 
necessarily actual) past, whether or not that past has 
been fully or partially realized. 

We have also shown at the end of the last section that 
in theories with a time reversal invariant Hamiltonian 
there are no time reversal invariant Lyapunov operators. 
In Ref. |9[ it was argued that in order to characterize a 
system as irreversible and single out a direction of time 
a Lyapunov functional should be time reversal invari- 
ant. Hence, if one accepts this view of Ref. Q then, 
by our result, quantum mechanics for finitely many par- 
ticles should indeed not be irreversible and should not 
exhibit an arrow of time if the Hamiltonian is time re- 
versal invariant. 



choice of generalized eigenvectors, {E\-4') = '0(^) is real. 
Then the first term on the r.h.s of Eq. (P5|) is the integral 
of a total derivative and therefore vanishes, as does the 
third term on the r.h.s. of Eq. ^, by Eq ^. Thus 

- (jdE\^\'^Y^hM^ . (Al) 

By Schwarz's inequality the last term can be estimated 
as 

i 

Y^ j dE\i:\'\h\' -Y^ j dE\^\''%\\ (A2) 



< 



where the first sum on the r.h.s. yields 1 and the equality 
sign holds if and only if 

b[{E)^-fh{E) , 7 = constant, (A3) 

which implies 

Since the 1. h. s. is purely imaginary, from Eq. ([M]) . this 
implies 7 = lA with A real. Thus, for real 4'{E), 



AT' > / dE\xP\E)\'^ 



(A5) 



with equality holding if and only if Eq. ()A3P holds with 
7 = iA , A real, i.e. if and only if 

b,{E) = Cic'^-^ , A real (A6) 

YME)ME) = ^|qP = 1. 
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Appendix A: Minimal variance and non-uniqueness 
of time operator 

^ We show for the case of a non-degenerate spectrum of 
H that minimal variance alone does not imply uniqueness 
of T. We first consider a state \ip) such that, with a given 



These functions give the same time operator and density 
as the single function 



biE) = e 



i\E 



(AT) 



With this choice AT' becomes minimal for real tp{E). 

For a state given by e^'^^-^^i/j{E), with real ipiE) and 
f{E), the same argument gives, upon replacing bi by 
^-iip{E) ^^^ i]xQi, one has minimal variance if and only if 



6,(^) = c,e'(^-^(^». 



(A8) 



This differs from Eq. dMl), as does the analog e*^^"'^^^)) 
of the single function in Eq. (jA7p . 

Hence among the set of all allowed functions bi{E) 
there is no choice of functions such that AT becomes 
minimal for all states with finite second moment. 
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